Time-dependent simulation of viscoelastic flows at high Weissenberg number using the log-conformation representation
نویسندگان
چکیده
We present a second-order finite-difference scheme for viscoelastic flows based on a recent reformulation of the constitutive laws as equations for the matrix logarithm of the conformation tensor. We present a simple analysis that clarifies how the passage to logarithmic variables remedies the high-Weissenberg numerical instability. As a stringent test, we simulate an Oldroyd-B fluid in a lid-driven cavity. The scheme is found to be stable at large values of the Weissenberg number. These results support our claim that the high Weissenberg numerical instability may be overcome by the use of logarithmic variables. Remaining issues are rather concerned with accuracy, which degrades with insufficient resolution.
منابع مشابه
Flow of viscoelastic fluids past a cylinder at high Weissenberg number: stabilized simulations using matrix logarithms
The log conformation representation proposed in [1] has been implemented in a fem context using the DEVSS/DG formulation for viscoelastic fluid flow. We present a stability analysis in 1D and attribute the high Weissenberg problem to the failure of the numerical scheme to balance exponential growth. A slightly different derivation of the log based evolution equation than in [1] is also presente...
متن کاملComputation of viscoelastic fluid flows at high Weissenberg number using continuation methods
The numerical simulation of viscoelastic fluid flow becomes more difficult as a physical parameter, the Weissenberg number, increases. Specifically, at a Weissenberg number larger than a critical value, the iterative nonlinear solver fails to converge, a phenomenon known as the High Weissenberg Number Problem. In this work we describe the application and implementation of continuation methods t...
متن کاملAn Adaptive Remeshing Strategy for Viscoelastic Fluid Flow Problems
The simulation of highly viscoelastic fluid flow problems is still plagued with the so-called high Weissenberg problem (HWP). The Weissenberg number We is the ratio of the magnitude of the elastic forces to that of the viscous forces. Most (if not all) numerical methods lose convergence i.e. fail to generate any result and/or fail to generate convergent solutions with respect to mesh size, at s...
متن کاملNumerical study of the square-root conformation tensor formulation for confined and free-surface viscoelastic fluid flows
We present a numerical study of a stabilization method for computing confined and free-surface flows of highly elastic viscoelastic fluids. In this approach, the constitutive equation based on the conformation tensor, which is used to define the viscoelastic model, is modified introducing an evolution equation for the square-root conformation tensor. Both confined and free-surface flows are con...
متن کاملA new numerical framework to simulate viscoelastic free-surface flows with the finite- volume method
A new method for the simulation of 2D viscoelastic flow is presented. Numerical stability is obtained by the logarithmic-conformation change of variable, and a fully-implicit pure-streamfunction flow formulation, without use of any artificial diffusion. As opposed to other simulation results, our calculations predict a hydrodynamic instability in the 4:1 contraction geometry at a Weissenberg nu...
متن کامل